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Abstract 

Stationary determinantal point processes are proved to be Brillinger mixing. 
This property is an important step towards asymptotic statistics for these 
processes. As an important example, a central limit theorem for a wide class 
of functionals of determinantal point processes is established. This result 
yields in particular the asymptotic normality of the estimator of the intensity 
of a stationary determinantal point process and of the kernel estimator of its 
pair correlation. 
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1 Introduction 

Determinantal point processes (DPPs) are models for repulsive point patterns, where 
nearby points of the process tend to repel each other. They have been introduced in 
their general form in [23] and extensively studied in Probability theory, see [HI and 
[23] • From a statistical perspective, DPPs have been applied in machine learning 
EDI. spatial statistics I221EI] and telecommunication PEI]. The growing interest 
for DPPs in the statistical community is due to their appealing properties: They 
can be quickly and perfectly simulated, parametric models can easily be constructed, 
their moments are known and the likelihood has a closed form expression. Their 
dehnition and some of their properties are recalled in Section ITTI and we refer to [22] 
for more details. Some realizations are showed in Figure [U 

We focus in this paper on stationary DPPs on the continuous space and we 
prove that they are Brillinger mixing. To the best of our knowledge, no mixing 
property was established so far for DPPs. The Brillinger mixing property is an 
important step towards asymptotic statistics for DPPs, which are mainly unexplored 
in the literature. The dehnition, recalled in Section 12.21 is based on the moments 
of the process. Specihcally, a stationary point process is Brillinger mixing if for any 
k >2 the total variation of its reduced factorial cumulant measure of order k is hnite, 
see for instance ]3] or [12]. Already known Brillinger mixing point processes include 
Poisson cluster processes and Matern hardcore point processes (of type I, type II and 
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some generalizations as in |32]), see na and [m. As far as we know, the Matern 
hardcore models are the only models of repnlsive stationary point processes that 
have been proved to be Brillinger mixing. Onr result shows that DPPs provide a 
new flexible class of repulsive Brillinger mixing point processes. 

In Sectional we give some applications of the Brillinger mixing property of DPPs. 
These are mainly based on general results established in PI, 0 and PI, that we 
extend and/or simplify in the setting of stationary DPPs. Namely, we prove the 
asymptotic normality of a wide class of functionals of order p of a DPP, in the spirit 
of m- This result allows in particular to retrieve the asymptotic behavior of the 
estimator of the intensity of a DPP, known since |2S], and to get the asymptotic 
normality of the kernel estimator of the pair correlation function of a DPP, which 
is a new result presented in Section 14.21 The Brillinger mixing property is useful for 
many other applications, see for instance m, nn and [IS]. In an ongoing project 
[2], this property is used to get the asymptotic normality of minimum contrast 
estimators for parametric DPPs. 

The reminder of this paper is organized as follows. Section [2] gathers some basic 
facts about stationary DPPs, moment measures of a point process and the Brillinger 
mixing property. Our main result stating that stationary DPPs are Brillinger mix¬ 
ing is presented in Section [31 Some statistical applications are given in Section jH 
Section El and Section El contain some technical proofs and Section [7] is an appendix 
dealing with the computation of the asymptotic variance in the statistical applica¬ 
tions of Section jH 

2 Preliminaries 

2.1 Determinantal point processes 

For d > 1, we denote by the class of bounded Borel sets on For x C 

and B G x(i?) stands for the number of points in x fl 5. We let J\f : = 

{x C M'^, x(i?) < oo, \/B G be the space of locally hnite conhgurations 

of points in R'^. This set is equipped with the cr-algebra generated by the sets 
{x C R'^, x{B) = n} for all B G Bo{W^) and all n G N U {0}, where N denotes the 
space of positive integers. A point process on R'^ is a measurable application from 
a probability space into the set Af. We denote a point process by a bold capital 
letter, usually X, and identify the mapping X and the associated random set of 
points. All considered point processes are assumed to be simple, i.e. two points of 
the process never coincide, almost surely. For further details on point processes, we 
refer to HE]. 

The factorial moment measures and especially the joint intensities of order A; of a 
point process, defined below, are important quantities of interest. They in particular 
characterize the law of determinantal point processes. 

Definition 2.1. The factorial moment measure of order k (k > 1) of a simple point 
process X is the measure on R'^^, denoted by such that for any family of subsets 
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Di,, Dk in 


(Di X ... X Dk) 


E 




^{xl€Dl,...,x^.eD^.} 




where E is the expectation over the distribution of X and the symbol 7 ^ over the sum 
means that we consider only mutually disjoints k-tuples of points Xi,... ,Xk. 

If admits a density with respect to the Lebesgue measure on this density 
is called the joint intensity of order k of IX. and is denoted by 

Important particular cases are the factorial moment measure of order one, called 
the intensity measure, and the factorial moment measure of order two. If X is 
stationary, then for all S C M'^, there exists p > 0 such that = PIS'!, where 

[S'! stands for the volume (Lebesgue measure) of S. In this case, for any x G M'^, 
(x) = p is called the intensity of the process and represents the expected number 
of points per unit volume. Regarding the joint intensity of order two, for (x, y) G 
and X y, p^‘^\x,y) may be viewed heuristically as the probability that there is 
a point of the process in a small neighbourhood around x and another point in a 
small neighbourhood around y. In spatial statistics, the second order properties of 
a point process are often studied through the pair correlation function (pcf). The 
pcf is dehned for almost every (x, y) G by 


9{x,y) 


p{x)p{y )■ 


In the stationary and isotropic case, p(x, y) = go{r) depends only on the Euclidean 
distance r = |x — p|. Intuitively, go{r) is the quotient of the probability that two 
points occur at distance r (taking into account the interaction induced by the pro¬ 
cess) and the same probability if there was no interaction. Consequently, for r > 0, 
a common interpretation, see for instance is that go{r) > 1 characterizes clus¬ 
tering at distance r while Po(^) < 1 characterizes repulsiveness at distance r. 

Determinantal point processes (DPPs) are dehned through their joint intensities. 
They have been introduced in their current form by Macchi in [23] to model the 
position of particles that repel each other. Since our results concern only stationary 
DPPs, we restrict the dehnition to this subclass, which simplihes the notation. 

Definition 2.2. Let C : ^ M. be a function. A point process X on is a 

stationary DPP with kernel C, in short X ~ DPP{C), if for all k > 1 its joint 
intensity of order k satisfies the relation 


p^^\xi, ...Xk) = det[C](xi,... ,Xfc) 

for almost every (xi,... ,Xk) G where [C](xi,... ,Xfc) denotes the matrix with 
entries C{xi — Xj), I < i, j < k. 

It is actually possible to consider complex-valued kernels and/or non-stationary 
DPPs, but this is not the setting of this paper and we refer to [TT] for a review 
on DPPs in the general case. The existence of a DPP requires several conditions 
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on the kernel C. Snfficient conditions in the stationary case are provided in the 
next proposition. They rely on the Fonrier transform of C and are easy to verify in 
practice, unlike the general conditions for non stationary DPPs, see [2]. 

We dehne the Fourier transform of a function h G as 


and extend this dehnition to by Plancherel’s theorem, see [SU]. We have the 

following existence result. 

Proposition 2.3 ([22]). Assume C is a symmetric continuous real-valued function 
in L'^(EA). Then DPP{C) exists if and only ifO< .F(C') < 1. 

In other words, by Proposition l2.3l anv continuous real-valued covariance function 
C in with < 1 dehnes a DPP. Henceforth, we assume the following 

condition. 

Condition K,{p). A kernel C is said to verify condition /C(p) if C is a symmetric 
continuous real-valued function in with C(0) = p and 0 < J^{C) < 1. 

By dehnition, all moments of a DPP are explicitly known. In particular, assum¬ 
ing /C(p), DPP{C) is stationary with intensity p and denoting g its pcf we have 


g{x,y) 


C{x-yf 

p2 


( 2 . 1 ) 


for almost every {x,y) G Consequently p < 1, which shows that DPPs exhibit 
repulsiveness. 

A hrst example of stationary DPP is the stationary Poisson process with intensity 
p, which corresponds to the kernel C{x) = pl{a;=o}- However, this example is very 
particular and represents in some sense the extreme case of a DPP without any 
interaction. In particular its kernel does not satisfy /C(p) since it is not continuous. 
In contrast, /C(p) is verihed by numerous covariance functions , and this makes easy 
the dehnition of parametric families of DPPs, where the condition J^(C') < 1 implies 
some restrictions on the parameter space. Some examples are given in [22] and [3], 
where the stationary Poisson process appears as a degenerated case. For instance, 
the Gaussian kernels correspond to C{x) = , x G M'^, where the existence 

condition implies a < l/(i/Fp^/'^). Another important example is the most repulsive 
stationary DPP with intensity p, as dehned and determined in [5]. Its kernel C is 
the Fourier transform of the indicator function of the Euclidean ball centered at the 
origin with volume p, which gives 


C{x) = 


^pF(f + 1) (2xAFF(f + l)^p^|x|) 


TT' 


d/4 


Vx G 


\X\ 2 


( 2 . 2 ) 


where Jd denotes the Bessel function of the hrst kind of order f. Some examples of 
2 ^ 

realisations of DPPs are given in Figure [H 
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Figure 1: From left to right, letting the intensity p = 100, realizations on [0,1]^ of a 
stationary Poisson process, a DPP with a Gaussian kernel and the maximal possible 
choice for the range parameter a {a = 0.056), a DPP with kernel (| 2 . 2 |) . 


2.2 Moment measures and Brillinger mixing 

In this section, we review the dehnition of the cumulant and factorial cumulant 
moment measures of a point process X as well as their reduced version. These are 
at the basis of the Brillinger mixing property dehned in the following. The relation 
with the Laplace and the probability generating functionals of X is also described. 
We assume that for any bounded set A, the random variable X( 74 ) has moments 
of any order. This ensures that the quantities introduced in this section are well 
dehned. Note that by dehnition this assumption holds true for a DPP. Further 
details on these topics may be found in [U E] and [H]. 

Definition 2.4. For k E N, the cumulant of the k random variables Xi,... ,Xk is, 
if it exists, 


Cum(Xi,...,Xfc) = 


dk 


dti... dtk 


logE 


exp UXi 


\i=l 


U=...=tfe=0 


The k-th order cumulant of the random variable X is Cumfc(X) := Cum(X ,... ,X). 
The notion of cumulant of random variables extends to point processes as follows. 


Definition 2.5. For k E N, the k-th order cumulant moment measure 7 ^ of a 
point process X is a locally finite signed measure on defined for any bounded 
measurable sets Ai,... ,Ak in by 

Ik fn = Cum 

/ \x€X xGX 



Definition 2.6. For k eN, the k-th order factorial cumulant moment measure 7 [fc] 
of a point process with factorial moment measure for r < k, is a locally finite 
signed measure on defined for any bounded measurable sets Ai,... ,Ak in R'^ by 

i[k] - !)• E { n 

^*=1 / \ki&Bi 
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where for all j < k, Vj denote the set of all partitions of {1,, k} into j non 
empty sets Bi,..., Bj. 

For stationary point processes, we can define the so-called reduced version of the 
previous measure. 

Definition 2.7. For any k >2, the reduced k-th order factorial cumulant moment 
measure of a stationary point process is a locally finite signed measure on 
defined for any bounded measurable sets Ai,... in by 

(n v)= fi Tfy (n (V - ^)) * 

where for i = 1,..., k — 1, Ai — x is the translation of the set Ai by x. 

The reduced cumulant moment measure is dehned similarly. An important prop¬ 
erty of signed measures is given by the following theorem leading to the dehnition 
of the total variation of a signed measure. 

Theorem 2.8 (Hahn-Jordan decomposition, see [H Theorem 5.6.1]). For any signed 
measure v, there exist two measures z/+ and z/“ uniguely determined by v such that 
at least one of them is finite and 


V = — V . 


Definition 2.9. Let u be a signed measure with Hahn-Jordan decomposition v = 
— v~. The total variation measure \ v\ of v is defined by 

\v\ = + v~. 


Following Theorem 12.81 for k >2, we denote the Hahn-Jordan decomposition of 
the reduced /c-th order moment factorial cumulant measure 

Definition 2.10. A point process is Brillinger mixing if, for k >2, we have 



< -l-cxo. 


The different moment measures of a point process X are related to the power 
series expansion of the Laplace and the probability generating functionals of X. 


Definition 2.11. The Laplace functional Lx of a point process X is defined for any 
bounded measurable function f that vanishes outside a bounded set of by 

Lx(/) = E(e-S^.x/(-)). 

Definition 2.12. The probability generating functional of a point process X is de¬ 
fined for any function h from into [0,1], such that 1 — h vanishes outside a 
bounded set, by 


G:s.{h) = E exp 


V log h{x) 




Notice that for any function h defined as in Definition 12.121 and taking values 
within a closed subset of ( 0 , 1 ], we have 


Gx(h) =Lx(-log(h)). 


Proposition 2.13 ([51 Section 9.5]). Let ^ be a point process with cumulant moment 
measures 7 ^ and factorial cumulant moment measures 7 [fc]. Let f and rj be bounded 
measurable functions on that vanish outside a bounded set. Assume further that 
rj takes values in [0,1]. Then, for all N E N, we have the following power series 
expansions when s > 0 and s —>• 0 


^ (-sV 

log Lx (s/) = 

j = l J' 

^ i-sV 

logGx(l-s7) = 51^-1— 
j=i J- 


J f{xi)f{x2 )... f{Xj)-fj{dXi 
J viXi)T]{x2) . ..T]{Xj)-f[p{dXi 


X dx2 X ... X dxj) + o(s^), 
X dx2 X ... X dxj) + o{s^). 


We conclude this section by giving the relation between 7 ^ and 7 [fc]. To this end, 
we recall the definition of the Stirling numbers of the first and second kind and refer 
to [H Section 5.2] for a detailed presentation. For a; G M and fc G M, we denote by 
x[k] — _ I'j _ _ l)l{o<fc<a;} the falling factorial of x. Assuming k < x and 

1 ^ J < ^5 the Stirling numbers of the first kind Dj ^ and of the second kind ^ 
are defined by the relations 


k k 

and x^ = y~] 
i=i i=i 

Proposition 2.14. Let A be a bounded set o/M'^. For any integer k, we have the 
relations 

7w(-4*) = 

3=^ 

k 

i=i 

Proof. We denote by L and G the Laplace and probability generating functionals of 
a point process with, for k > 1, factorial cumulant moment and cumulant moment 
measures 7 [fc] and 7 ^, respectively. By Proposition 12.131 for all A G N, we have as 
s —>• 0 , s > 0 


N ( \k 

logG'(l - s 1 {, 6 A}) = [E 

k=i 


(2.3) 


As noticed after Definition 12.121 

logG(l - sljseA}) = logL(-log(l - sljseA})) = logL(-log(l - s)l{sg 7 i}). 
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Since s 


r\j 


— log(l — s) as s —)■ 0, we have by Proposition 12.IHL 


logG(l - sl,,e.4)) + ois"). 

1=1 2- 


By [H (24.1.3.1.B)] we deduce that 



N { \k k 

= E ^ + o(s"). (2.4) 

k=l j = l 


We conclude by identifying the coefficients in fl2.3p and 02.41) . The proof of the second 


formula is similar, starting with the other powers expansion in Proposition 12.131 and 
using P, (24.1.4.1.B)] instead of P (24.1.3.1.B)] 


□ 


3 Main result 

In this section, we prove in Theorem 13.21 below that a DPP with kernel verifying the 
condition IC{p) is Brillinger mixing. We recall that this mixing property involves the 
factorial cumulant moments of the DPP. It is not easy to deduce these moments from 
the initial Dehnition 12.bl However, the power series expansion of the log-Laplace 
functional in Proposition 12.131 which is known for a DPP, allows us to derive a 
closed form expression for the factorial cumulant measures as stated in the following 
lemma. 

Lemma 3.1. Consider a DPP with kernel C verifying eondition /C(p) and, for k G 
N, denote its k-th factorial cumulant moment measure by 7 [a;] . For every measurable 
bounded set A in and k >2, we have 


'^[k\{A^) = (-l)^+^(/c - 1)! [ C{x2 - xi)... C{xi - Xk)dxi ...dxk 


Proof. By [2B1 Proposition 3.9], we deduce that for any bounded set H C and s 
small enough. 






Then, by Proposition 12.131 we have by the last equation that for all p G N and any 
bounded set H C 


p 


p\ 




n ■pi!p2! ■ ■ -Pn 


n=l 


Pl + ..-+Pn=P 
P1,...,P„>1 
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Thus, by Proposition 12.141 we have for k>2, 


7w(^‘)=i:(-i) 

p=l 


k— 




pi 


n=l 


ri- Pi!p2!---Pn! 


/ C{x 2 — Xi) . . .C{xi — Xn)dxi . . . dXn- (3.1) 

JA^ 


By [H (24.1.2.I.B)], it is easily seen that 


E 


pi 


= E 


pi 


E 


pi 


■ • -Pn'- p,+.^p„=pPl'-P2'- ■ ■ -Pn'- p,+...r;^_,=p PlW- ' ' ' Pn-P 

Pl,..;Pn>l 

= nP-{n- 1)P (3.2) 


By dehnition 


J2{-l)^^-PDp^kinP - (n - 1)P) = - {n - 1)1^1 

p=i 


(3.3) 


which is null for every n < k. Therefore, by fl3.2p and 03.31) . only the terms n = k is 
non null in the sum 03.11) . □ 

We are now in position to prove our main result. 

Theorem 3.2. A DPP with kernel verifying the condition K,{p), for a given p > 0, 
is Brillinger mixing. 

Proof. For any t > 0, we have by taking / = in Definition 12.71 

7[fc]([-t,t]'^^) = dx. 

By Lemma 13.11 

y dx = (-l)^+^(fc - 1)! hit) (3.4) 


where for all fc > 1 and f > 0, Ik(t) := C{x 2 — Xi ).. .C{xi — Xk)dxi... dxk- 

Since C verifies the condition /C(p), by Mercer’s theorem, see also PH Section 2.3], 
we have for alH > 0, 

C{x-y) = '^Xj{t)(j)j{x)(j)j{y), W{x,y) e [-t,t]‘^, 
jm 

where for all j G N, {0j}jgn is an orthonormal basis of and Xj{t) belongs 

to [0,1] by [231 Theorem 4.5.5]. Then, by orthogonality of the basis we have 

for alH > 0 and k > 1, 


4(t) = E\'W<E\W = ^iW (3-5) 

jgn jeN 
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where Ii{t) = J^_^ ^dpdx = 0{t'^). Thus, by Theorem I2.8[ fi:h4l] and fIThp . there 
exists a constant k > 0 and T > 0 such that for all t >T, 










k-1 


dx 


< Kt^. 
(3.6) 


Henceforth, we assume t > T. By Theorem 12.81 at least one of the measure 
or is hnite. Let us assume without loss of generality that is hnite. Thus, 
by fid.hi) and the monotonicity of the measure we have 

P, ((!“*■ ‘I' ■ + 27,^7((K‘‘)‘-‘)) , (3.7) 

so by positivity of 

L.r dT(('-*’‘1'-d‘^‘)2((!-*.*1'-d‘^‘)(3.8) 


Further, for all {x, y) G 


—t t 
2 ’ 2 


1 d 


nL L 
2 ’ 2 


n2d 


, y + X ^ SO for all x G 


nL L 
2 ’ 2 

-\-red 


1 d 


we have 


C [—t,tf — X. It follows by fl3.8p and the monotonicity of that 


:7[fc] 


+red 


i-t 

Y'2 




Hence by ([221)-dSS]), we have 

-t 


+red 

7[fc] 


2 ’ 2 


< 


(k + 2‘'7[-ppM''('=-d)) 


By letting t tend to inhnity in the last equation, we see that is hnite and so is 


red 


7[fc] 


by Dehnition 12.91 which concludes the proof. 


□ 


4 Statistical applications 

Many applications of the Brillinger mixing property for point processes may be 
found in m, CO], m, mi and Ba. We present in this section some of these 
applications for DBFs. We prove in Section 14.11 a general central limit theorem 
for certain functionals of a DPP that are involved in the asymptotic properties of 
standard estimators. As an example, we apply this result to the estimator of the 
intensity of a DPP. Another important application concerns the asymptotic behavior 
of minimum contrast estimators for parametric DPPs, which will be the subject of 
a separate paper. In Section 14.21 we obtain the asymptotic properties of the kernel 
estimator of the pcf of a DPP. In particular, we prove a central limit theorem for 
the pointwise estimator of the pcf and for its integrated squared error. 


10 























4.1 Asymptotic behaviour of functionals of order p 

We present an important consequence of the Brillinger mixing property, namely a 
central limit theorem for a wide class of functionals of the point process and the 
convergence of their moments. A hrst theorem was mentioned in na and proved 
in [ig. We present here a more general version that yields in particular the asymp¬ 
totic normality of standard statistics as the natural estimator of the intensity of the 
process. These results apply to stationary DPPs under condition /C(p) as explained 
and exemplihed at the end of this section. 

For a given set D of we denote by dD the boundary of D. 

Definition 4.1. A sequence of subsets {D„}neN of is called regular if for all 
n eN, Dn C Dn+i, Dn is compact, convex and there exist constants ai and such 
that 

ain'^ < \Dn\ < 

< Hd-i (dDn) < a2n'^~^ 

where T-Ld-i is the {d — 1)-dimensional Hausdorff measure. 

Note that any sequence of subsets as above grows to in all directions. For 
p > 1, let /ij be a function from into R that depends on a given set D C R*^ 
and dehne for a stationary point process X, 

Np{fD)-= Y. fD{Xl,...,Xp). 

{xi,...,Xp)£X.P 

By letting the set D in the last equation be a sequence of regular subsets 
we have under some suitable conditions on the function fo^, the following central 
limit theorem on the sequence {Ap {fD„)}n£w proof is postponed to Section[S21 

Proposition 4.2. Let and {Dn}ne¥i be two sequences of regular sets in the 

sense of Definition \4.1\ such that k for a given n > 0. Assume that there 

exists a bounded and compactly supported function F from into R"*" such that 

for alln eN and [xi,... ,Xp) E R'^^, 

\fDSxi,---,Xp)\ < Xi,...,Xp-Xi). (4.1) 

\^n\ 

Assume further that the point process X is ergodic, admits moment of any order and 
is Brillinger mixing in the sense of Definition \2.lU . Then, for all k>2, we have 

Cum^ y\I\\N, UdS) = O . (4.2) 

Moreover, if there exists a > 0 such that 

Var (dFlA'p (/dJ) (4.3) 

we have the convergence 

mO,X) ( 4 , 4 ) 

and the convergence of all moments to the corresponding moments o/jV(0,(T^). 
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By fl4.2l] . the variance given in fl4.H|) is uniformly bounded with respect to n G N. 
If Dn and in fl4.3p are sufficiently generic, the convergence 04.31) of the variance 
holds true. However, in the general case, it must be assumed. To check 04.3p 
in applications, it is convenient to express the variance in 04. 3 p in terms of the 
factorial cumulant moment measures of X. In appendix, we detail this expression 
for the important situations p = 1 and p = 2 with 2 : 2 ) = 0 for xi 7 ^ X 2 , see 

Lemmas 17.1117.21 and 17.31 

Proposition 14.21 applies to stationary DPPs with kernel verifying /C(p) provided 
04.11) is verihed. Indeed, Soshnikov in [2B] proved that a stationary DPP is ergodic. 
Moreover, a DPP admits moments of any order by dehnition and is Brillinger mixing 
under condition /C(p) by Theorem 13.21 As a direct application when p = 1, we 
retrieve a result of giving the asymptotic normality of the estimator of the 
intensity of a DPP. 


Corollary 4.3. Let ^ be a DPP with kernel verifying IC{p) for a given p > 0 and 
be a family of regular sets. Define for all n G N, 


Pn 


I n I ^ l{a:6nin}- 
xex 


(4.5) 


We have the convergence 



p) ^^iV( 0 ,a 2 ) 

n—>-+oo 


where = \\mn^+ooV ar (^^\Dn\pn) = P - fud C(x)^dx. 

The proof of this corollary follows by taking p = 1 and foS^) = \^\^{x£Dn} 
in Proposition 14.21 In this case, the assumption 04.31) holds by Lemma 17.11 and a 
straightforward calculus. 


4.2 Applications to the empirical pair correlation function. 


We consider in this section the estimation of the pcf of a stationary and isotropic 
DPP in M'^. In this setting g{x,y) = go{r) depends only on the Euclidean distance 
r = \x — y\. Let be a sequence of regular subsets of in the sense of 

Dehnition 14.11 { 6 n}neN a sequence of positive real numbers, and k a function from 
M into For n G N and z G M'^, we denote for short Df := Dn — z the translation 
of Dn by z. For r > 0, we consider the kernel estimator of go{r) 


9n{r) 




l{xen>n, 

{x,y)&yP 

x^y 


y^Dn} 


1 . / r - [a; - p| \ 

bnlDnHDr^l [ bn ) 


(4.6) 


where pn is given by fl4.5|) and cr^ = denotes the surface-area of the d- 

dimensional unit sphere. Some comments and details about this estimator may 
be found, for instance, in [2S1 Section 4.3.5] or [S]. 

The following proposition gives the asymptotic normality of the pointwise esti¬ 
mator p„(r) for r > 0. Its proof, given in Sectional is based on Proposition 14.21 and 
results from HD]. 
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Proposition 4.4. Let be a regular sequence of subsets of Assume 

that the sequence {6ri}neN such that b^\Dn\ —)■ +cx3 and &^|-Dn| 0. Let k be 

a symmetric and bounded function with compact support included in [—T, T], for a 
given T > 0, and J^k{t)dt = 1. Let C be an isotropic twice differentiable kernel 
on \ {0} verifying IC{p) for a given p > 0. Then, for all r > 0, we have the 
convergence 


^/bn\Dn\ idnir) - go{r)) iV(0, r^) 

' n —>-+00 

where rf: = J^k^{t)dt. 


In addition to the previous result, we state the asymptotic normality of the 
integrated squared error of the estimator pfign where Pn is dehned in fl4.6p . This 
quantity is the basis of an asymptotic goodness-of-£t test for stationary DPPs as 
presented in [5]. For all segment I C \ {0} and n G N, denote 

ISE„(J) = ^ {pn 9 n{r) - p‘^go{r)y dr. 

Proposition 4.5. Let {-DnjneN be a regular sequence of subsets of Assume that 
the sequence {6n}neN is such that 6^ —t 0 and bn\Dn\ —)■ +cxo. Let k be a symmetric 
and bounded function with compact support included in [—T, T], for a given T > 0, 
and J^k{t)dt = 1. Let C be an isotropic twice differentiable kernel on \ {0} 
verifying JC{p) for a given p > 0. Then, for all segment I C M’*' \ {0}, we have as n 
tends to infinity, 

6„|DjE(ISE„(J)) = 2p2 f -^^dr [ k{tfdt + Oibr,)+0{\Dfibl). 

Ji (Jdr^~ Jr 

If in addition h^f\Dn\ —)■ 0 then 


^n\Dn\ (ISE„(J) - E (ISE„(J))) iV(0, r^) 

n. —^4-00 


where 

uct. 



(k * kfi(s)ds and * denotes the convolution prod- 


Proposition 14.51 is an application to the DPP’s case of the results given in [S]. 
In addition to the Brillinger mixing, ensured by Theorem ld.2[ and the properties of 
the sequence the authors need two additional assumptions. Namely, these 

assumptions are the locally uniform Lipschitz continuity of the first derivative of 5^0 
and a second assumption related to the densities of the reduced factorial cumulant 
measures. By fl2.ip and since C is twice differentiable on M'^\ {0}, the hrst derivative 
of go is uniformly Lipschitz continuous on every compact sets in R"*" \ {0} so the 
hrst assumption holds. The second assumption is verihed by Lemma 14.bl below. 
Consequently, Proposition 14.51 is proved by [3 Lemma 3.4] and [HI Theorem 3.5]. 
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Lemma 4.6. Let be an isotropic DPP with kernel C verifying the condition K,{p), 
whose reduced factorial cumulant moment measures of order 3 and 4 have densities 
and c|’ 4 p respectively. For all compact set K and e > 0, we have 


sup 

(|«|>|) g ( K ®')2 


red, 
C[3] I 


u,v) 


< +00 


(4.7) 


and 


sup 

(| n |>|)€( X ®^)2 


Cuf{u, w,v + w) dw < +CX), 


14] 


(4.8) 


where = K + 5(0, e) and 5(0, e) is the Euclidean ball centred at 0 with radius 

e. 


Proof. By fl7.2^ - 07.31) in Section [71 we have for all {u,v,w) G 

c\f({u,v) = 2C{u)C{v)C{v - u) 

and 


Cyff{u, V, w) = —2 [C{u)C{v)C{u — w)C{v — w) 

+C{u)C{w)C{u - v)C{v -w)+ C{v)C{w)C{u - v)C{u - w)]. 

Notice that K®'^ is compact and since C verihes the condition /C(p), it is continuous. 
Therefore, by 07.21) . 04.71) holds immediately. Finally, 04.81) is verihed by Cauchy- 
Schwarz inequality and 07.31) . □ 


5 Proof of Proposition 14.2 

5.1 Complement on the moments and cumulants of a point 
process 

We present here the necessary background to prove Proposition 14.21 Let p and k be 
two integers and X a point process that admits moments of any order. Consider, 
for 1 <i <k, the random variables 

Np{(t)i)= Y. (i)i{xi,...,Xp) ( 5 . 1 ) 

{xi,...,Xp)&X.P 

where for i = 1 ,... ,p, 0, is a function from to M. 

For l,s < kp , denote (resp. Q[) the set of all partitions of {l,...,kp} 
(resp. {1,..., /}) into I (resp. s) non empty sets pi,... ,pi (resp. qi,..., qs). For 
r = 1 ,..., s, denote 01 ,..., (3\qp the elements of the set qr and |g| the cardinal of a 
given set q. Then, as proved by Jolivet in [ini pl21-122], we have 

kp I 

= Y. T. H ( 5 . 2 ) 

^=1 n,G(Pf^ xieQl 
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where for all /, s < kp, 


/im.xi) 



k s 

X X{ct)i{d(^i-i)p+i,... ,dip) (5.3) 

i=l r=l 


The introduction of the term 6 is not easy to understand at hrst sight. For the sake 
of clarity, we give an example for p = k = 2 and 112 := {pi,P 2 } a given partition 
of the set {1, 2, 3,4} into 2 non empty sets, namely pi = {1,4} and p 2 = {2, 3}. In 
this case, we have 


k 

n 0 i(^(*-l)p+l; ■ ■ ■ 5 d2)4>2{d3, 6 ^ 4 ). 

i=l 

Thus, by the last equation, we have 
2 k 

n n ^{^m=dj} n i)p+i! ■ ■ • ^ “ 01 ( 2 ^ 1 ) 2 ^ 2 ) 02 ( 2 ^ 252 ^ 1 ) 

m=lj£pm i=l 

and a similar calculus is done if, for / = 1,.. ., 4, we choose another partition If; of 
{1, 2, 3, 4}. We can now describe completely Cum{Np (0i ),... ,Np (0^)). 

Theorem 5.1 ([IS])- The cumulant moment Cum(iVp (0i ),,Np (0^)) is equal to 
the sum of integrals in Formula fl5.2|) that are indecomposable, i.e. that 

can not he decomposed as a product of at least two integrals. 


5.2 Proof of Proposition 14.2 

Assuming (14.21) and (14.31) . the proposition is proved by [ISl Theorem 1 ]. Let us check 
(14.21) . By [271 Chapter II, Section 12, Equation (37)], if X and Y are two independent 
random variables Cumfc(X + T) = Cumfc(X) + Cumfc(y) and the cumulant of order 
A; of a constant is null for k > 2. Consequently, for k >2, 

Cum, [Np {fnj - E {Np (AJ)]) = Cum, A (Aj) 

= |ApCum,(A(/Dj). 


By Theorem 15.11 for every A: G N, Cum, {Np (/o^)) is a hnite sum of indecom¬ 
posable integrals Ii{Ili,xi)- Thus, it is sufficient to prove that for any k > 2, each 


integral A(A,xi) 


= 0 \D 


\l-k 


By (15.31) we have 


„ l k 

7071/, y;^) — I ^ n n 7{a;„i=6»j} /D„(^(j-l)p+l5 ^(i-l)p-|-25 • • • 5 ^ip) 

® m=ljepm i=l 


X n A( 72 :/ 3 i... (5.4) 


r=l 
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Then, by Definition I2.9[ we obtain from the last eqnation that 


I k 

\h(n.l,xi)\ < [ n n '^{Xm=9,}Yl\fDr,{0{i-l)p+l,0{i-l)p+2, ■ ■ ■ ,0ip) 
® m=ljGpm * = 1 




r=l 


Using fl4.ip . we get 


i-^;(nnx<j)i —1~ 1^ rr n ^{^m=6j}rTi{0{i_i)p+ie-Dn} 

|i^n| m=ljepm i=l 

s 

X F 0(j_i)p+i, • • •, dip \lqr \ {dxp^... (5.5) 

r=l 

Let ||U||oo denotes the snpremnm of F on Since the fnnction F is 

bonnded and compactly snpported, there exist compacts iLi,..., Kp_i snch that 

V(a7i, . . ., Xp—i) G (M , F{xi^ ■ ■ ■ 19yp—i) ^ I IT’I . .. l{xp_igiCp_i}- 

Then, we dednce from flh.hp that 





k 

} n ^{9(i-l)-p+l&Dn} 
2=1 


p-1 

IT —l)p+?? + l 

p=l 


— ®(i-l)p+l)SiCT7} n \lqM.dxp^...dxfi^ 


r=l 


Qr I ^ 


(5.6) 


Moreover, as already proved in [13 Section 4, Theorem 3], we have as n tends to 
infinity. 



k 

} IT ^{^(i-l)p+l€n>n} • • • 

2=1 


P-1 

■ IT ^{{^(i-l)p+r;+l 

? 7=1 


n \lqr\{dx0^... = O {\Dn\). (5.7) 

r=l 


Since - > k, the right hand term of fl5.6p is, by 05.71) . asymptotically of order 

n^ + OO 

\Dn\^~^, which ends the proof. 


6 Proof of Proposition 14.4 

The proof is based on the following lemmas. 

Lemma 6.1. Let {Dn}nm be a regular sequence of subsets of M'^. Assume that the 
sequence {&ri}neN is such that —)• 0 and —t +oo. Let k be a symmetric and 

bounded function with compact support included in [—T, T], for a given T > 0, and 
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/]g k{t)dt = 1. Let C he an isotropic twice differentiable kernel on \ {0} verifying 
/C(p) for a given p > 0. Then, for all r > 0, we have the convergence 


\/bn\Dn\ {pn9n{r) - E(p2^„(r))) 


n—>- + 00 


where f^k‘^{t)dt. 

Lemma 6.2. Under the same assumptions as in Proposition for all segment 
I C R’*' \ {0}, there exists a constant M > 0 such that 


sup 

re/ 


E {Plgnir) 


P^ 9 o{r)) <blMp^ f t^\k{f)\dt. 

«/M 


The proofs of Lemmas I6.HI6.2I are postponed to the end of this section. Let us 
now prove Proposition 14.4[ For all n G M and r > 0, we have 


Pn 


- go{r)) = An + Bn + C„ 


where 


-^n. \/bn\D' 


Bri 


Pndnir) -E{p^gnir)) 

Hpn9n{r)) - p'^goir) 


Cn = Jbn\Dn\goir) p^ - p^ 


By Lemma 16.11 we have the convergence 


A distr. T^T/r\ 2 

An -)■ A^( 0 , K 

n—>-+cxD 


( 6 . 1 ) 


( 6 . 2 ) 


and since bn\Dn\ tends to 0 as n tends to inhnity, we have by Lemma ESI 


Bn 0. (6.3) 

n—>-H -cxd 

By Corollary 14.31 and the delta method, we know that ^J\Dn\{Pn — P^) converges in 
distribution. Since 6^ —)■ 0, we deduce that 

Cn 0. (6.4) 

n—>-+oo 

Finally, by inserting fl6.2p - fl6.4l) in fl6.ll) . the proposition is proved by Slutsky’s the¬ 
orem and the almost sure convergence of p^ to p^. 


6.1 Proof of Lemma 16.1 

We need the following result. 

Lemma 6.3. Let r > 0 and D a subset of such that the Euclidean ball B{0,r) 
is included in D. Then, for all x G 5(0, r), we have 5®^ C 5 fl . 
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Proof of Lemma 1 6'., Vi By definition, D fl = {u E D, u + x E D} and 

= {u E D, \/v E B{0, r), u + v E D} . 

Since x E B{0,r), we have the inclnsion C -D fl D^. 

Dehne for all n G N and {xi,X 2 ) E 


□ 


fDr,iXl,X2) = 1 {xiGDn, X2eDn} 




rk 


' r — \xi — X2I 


Notice by 04.61) that 


inffjr'' V3n(f) = Y1 fDn(Xl,X2). 


(6.5) 


{xi,X2)€X^ 

Xl^X2 


The snpport of k is inclnded in [—T, T] so for any (xi,X2) E 


t)2d 


r — Xi — X2 


l{x2Gn>„} < 


< 


k 


r — Xi — X2 


r — Xi — X2 


^{\x\-X2\<r+Tbn} 


^{\x\—X2\<r+T} 


as soon as 6^ < 1 which we assnnie withont loss of generality since {6n}nGN tends to 
0. Then, by Lemma [6.31 and since k is bonnded, there exists M > 0 snch that for 
all n G M, 

ifonixii X2)\ < Q(^r+T): ^{\xi-X 2 \<r+T}- 

\ \ 

Therefore, by 04.2p in Proposition 14.21 and 06.Sp . we have for all k>3, 

Cnmfc (^^\Dn\bn(Tdr'^~^Pn9n{r)^ = 0{\Dn\^~^) 
whereby for all A; > 3, 

Cnmfc (^^Jbn\Dn\p^gnir)^ = O (b~^^'^\Dn\^~^^ 

which tends to 0 when n goes to inhnity since b'^\Dn\ -E 00 . Further, the conver¬ 
gences of for /c = 1,2 are proved in (TU] under conditions 

that we have already verihed after Proposition 14.51 Finally, Lemma [6.II is proved by 
the cumulant method, see m Theorem 1]. 


6.2 Proof of Lemma 16.2 


By 04.6p and Dehntion 12.11 we have for all n G N and r G /, 

.2 r r 1 


E {Pn9n{r) 


p 


^{x£Dn,y€Dr, 


crdr^~^bn 


p 


D„ n 


^{x£Dn,yGDr, 


adV^-^br, 


Dr, n D 


x-y 


\ br, J p^ 

k (- -^ 9o{\x - y\)dxdy. 
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To shorten, denote kb„{.) = By the snbstitution 2; = x — y and since y ^ 

if and only if z G D^, we obtain from the last eqnation that 


E {pndnir)) 




CTdr 

P^ 


d-l 


\Dn^Dl\ 


h^r - \z\)go{\z\)dzdy 


(JdT 


d-l 


hAr - \z\)go{\z\)dz. 


Converting this integral into polar coordinates and by symmetry of k, we get 

n+oo <^-1 


E {PnMr) 


= p 


= p 


hAr - t)go{t)dt 




d-l 


go{r + ubn)du. 


For all n large enongh, we have for all r G / that 7^ > T, hence 

On 


E {pl9n{r)) = j^^k {u) 


V + ubrt 


d-l 


go{r + ubn)du. 


Assnme that I writes [rmin-, rmax] for Vmax > ^min > 0 and dehne for s G M’*', 
/(s) := go{s). Notice that := [r^m - Tbn,rmax + Tbn] C M+ \ {0} as 

soon as n is large enongh which we assnme withont loss of generality. Since 5'o(-) is 
of class on so is /(.). Thus by Taylor-Lagrange expansion, we have 


^{pl9n{r)^ = P^ j ^k{u){^f{r) + f\r)ubn + J f"{s){ubn + r - s)ds'j du. 

( 6 . 6 ) 

Since k is symmetric, we have j]^j-uk(u)du = 0. Moreover, 


sup |/"(s)| < sup 


d-l 

min 




showing that f"{.) is uniformly bounded on by a constant M. Further, for 

all n G N and s G [r,r + ubn], \ubn + r — s| < \ubn\- Finally, since j]^j,k{u)du = 1, 
by 06.60 . we obtain for n large enough, 

E (Pndnir) - P^9oir)) < blMp^ f u^\kiu)\du, Vr G I. 

«/M 


7 Appendix 

We gather here some results useful to compute the asymptotic variance in Proposi¬ 
tion Upland Corollary 14.31 

Let X a stationary point process on and c^p and the densities of 
its factorial cumulant moment measures of order 2, 3 and 4, respectively, assuming 
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they exist. If X is a DPP with kernel C verifying the condition IC{p), for a given 
p > 0, then we deduce from Dehnitions 12.21 and 12.61 that for all {u,v,w) G 



-C\u), 

2 C{u)C{v)C{v-u), 


- 2 


C{u)C{v)C{u — w)C{v — tc) + C{u)C{w)C{u — v)C{v 


+ C{v)C{w)C{u - v)C{u - w) 


(7.1) 

(7.2) 
w) 

(7.3) 


Lemma 7.1. Let f be a function from into M that is bounded, measurable and 
compactly supported. Then we have 


Var ^ /(x) 




f{x)f{x + y)c\f^^{y)dxdy + p [ f{x)dx. 


Proof. Notice that 


Vxex / (a:,y)ex2 xex 

Then, denoting the density of the second order factorial moment measure, we 
have by Dehnitions 12.11 and 12.hi 


Var ( ^/(x) ) = / f{x)f{y) (p^‘^\x,y) - p'^) dxdy + p [ f{x)dx 
= f f{x)f{y)ci 2 ]{x,y)dxdy + p f f{x)dx. 

7M2d J^d 

Finally, by Dehnition 12.71 we have 

Var (^ f{x)] = f f{x)f{x + y)c\ff{y)dxdy + p f f{x)dx. 
VxGX / 


□ 

Lemma 7.2. Let f be a function from into M that is bounded, measurable and 
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compactly supported. Then, we have 


Var I ^ f{x,y) 

y{x,y)eX.^ / 

= {fix, x + y) + fix, X + y)fix + y, x)) c\^fiy)dxdy 
+ j {fix, y) + fx, y)fiy, x)) dxdy 

JR2d V / 

+ / [fix, x + y) + fix + y, a:)] [fix + y,x + u) + fix + u,x + y)] c\fiy, u)dxdydu 

JR3d i j 

+ 2p / [fix, y) + fiy, a:)] [fiy, y + u) + fiy + u, y)] c\f^iu)dxdydu 
+ P [fix, y) + fiy, x)] [fiy, x + u) + fix + u, y)] c\fiu)dxdydu 

JRSd ^ J 

+ P^ f, [fi^^ y) + fiy^ 2;)] [fiy, u) + fiu, y)] dxdydu 

JR3d 

+ / fix,x + y)fix + u,x + v)c\fiy,u,v)dxdydudv 
7R4d ^ J 

+ 4p / fix, y)fiy + u,y + v)c\fiu, v)dxdydudv 

aR4d ^ ^ 

+ 2 J f{x, y)fix + u,y + v 
+ 4p^ / fix,y)fix + u,v)c\fiu)dxdydudv. 

aR4d ^ ' 

Proof. This lemma is a generalization of PEI Lemma 5] for a function / non necessary 
symmetric. The variance is first computed with respect to the factorial moment 
measure by Definition 12.11 Then, the factorial moment measure is written in terms 
of the factorial cumulant moment measure by [H Corollary 5.2 VII] and the result 
is obtained by using Definition 12.71 We refer to the proof of PEI Lemma 5] for the 
detailed calculus, the only change being the use of the following decomposition in 
place of the original one, 

7^ 


I] fix,y)] = ff:y) + ff:y)fiy:f 


d^.?/)ex2 


(a:,y)eX2 

+ H [fix,y) + fiy,x)][fiy,u) + fiu,y)] 

(x,y,w)GX^ 

+ ff,y)fiu,v). 


□ 


Lemma 7.3. Let f be a function from into M that is bounded, measurable, 
symmetric and compactly supported. Let h be a function from into M that is 
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hounded, measurable and compactly supported. Then, we have 


Cov I ^ f{x,y),'^h{u] \ = 

K.{x,y)GX^ «GX 


f{x, X + y)h{u + x)(f^^{y, u)dxdydu 


+ p f{x, y)h{u + x)c\2]{u)dxdydu 

JRSd I 1 

+ p f{x, y)h{u + y)c\2]{u)dxdydu 

JRSd I J 

+ / f{x, y + x) [h{x) + h{y + x)] c\^f{y)dxdy 

+ y) [Kx) + h{y)] dxdy. 


Proof. Notice that 




H /(a;>2/) I ^(«)) = H f{x,y)h{u)+ Y. f{^^y){Kx)+ h{y)). 

,y)eX2 / \«GX / (x,y,u)GX3 (x,y)€X^ 


Then, by the last equation and Definition 12.11 we have 
Cov I ^ f(x,y),Y I = 

\(x,y)GX^ ueX 


f{x, y)h{u) (p® (x, y, u) - pp^‘^'> {x, p)) dxdydu 

JJ fix,y)ih{x) + h{y))p^^\x,y)dxdy. 


Finally, the proof is concluded by in Corollary 5.2 VII] and Definition 12.71 


□ 
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